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Abstract. A two-dimensional lattice gas model is proposed. The gromid state of this 
model with a fixed density is neither periodic nor quasi-periodic. It also depends on 
system size in an irregular manner. On the other hand, it is ordered in the sense 
that the entropy density is zero in the thermodynamic limit. The existence of a 
thermodynamic transition associated with such irregularly ordered ground states is 
conjectured from a duality relation for a thermodynamic function. This conjecture is 
supported by a phcnomcnological argument and numerical experiments. 



PACS numbers: 64.60.De,05.60.-hq,75.10.HK 
1. Introduction 

The ground state of a crystal is described by the repetition of a unit cell. Recalling 
that the repetition is a simple deterministic rule generating an infinite size pattern 
from a finite size pattern, we become curious about equilibrium systems with ground 
states described by other deterministic rules than the repetition. Note that the entropy 
density of such deterministic ground states is zero in the thermodynamic limit. In this 
sense, the ground states are ordered as observed in crystals. We call them irregularly 
ordered ground states. In this paper, we are concerned with thermodynamic transitions 
associated with irregularly ordered ground states. 

We shall give a precise definition of irregularly ordered ground states by considering 
some examples. First, quasi- crystals provide aperiodic ground states [1]. However, since 
the Fourier transform of a quasi-crystal configuration consists of sharp delta peaks, we 
classify them as regular ground states, following the terminology in studies of dynamical 
systems [21 [3]. Another type is represented by ground states of the anti-ferromagnetic 
Ising model on a triangular lattice, which are often called disordered ones. The ground 
states are neither periodic nor quasi-periodic. However, since their entropy density is 
not zero as calculated explicitly we do not regard such disordered ground states as 
irregularly ordered ground states. It is worthwhile noting that a string representation 
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of the ground states is described by a stochastic rule corresponding to non-intersecting 
random walks [5l[6]. Based on these considerations, we define irregularly ordered ground 
states as those which have no delta peaks in their Fourier transform and whose entropy 
density is zero. 

As an example with irregularly ordered ground states, one remembers the Frenkel- 
Kontorova model O [8], because its ground state is described by a chaotic map. Since 
it does not exhibit a phase transition at finite temperature, it is natural to consider an 
extended Frenkel-Kontorova model in two- dimensions. Recently, for such an extended 
model together with other related models, a certain type of transition has been found 
with the construction of ground states [9]. An interesting study with experimental 
relevance is now going on. However, from a theoretical viewpoint, the model is still 
complicated if one attempts to elucidate the nature of thermodynamic transitions 
associated with irregularly ordered ground states. In order to develop a theoretical 
argument, it seems better to study a two-dimensional lattice model. We here remark 
that non-trivial ground states in lattice gas models, which include an infinitely series 
of periodic ground states, were discussed in Refs. [TOl HI]. Furthermore, we should 
mention that Ref. [I2| reports a three-dimensional short-range interaction model in 
which non-periodic long-range order sustains at positive temperature. This reference 
also provides a summary of mathematical results on non-periodic order of ground states 
and Gibbs states for lattice models. 

Based on these previous studies, in this paper, we first present a two-dimensional 
lattice gas model with irregularly ordered ground states. The ground states of our model 
are described by a one-dimensional cellular automaton whose deterministic rule yields a 
non-repetitive structure given by a superposition of Sierpinski gaskets. We next consider 
a possibility of thermodynamic transitions associated with the ground states. Since our 
model turns out to be self-dual for a duality transformation, we can present a simple 
argument for the thermodynamic transition. 



We consider hard-core particles in a triangle lattice consisting of sites. Let A be the 
set of all sites in the lattice and A* be the set of all upward triangles. We express a 
configuration of the particles by an occupation variable cxj for i G A: = 1 when site i 
is occupied by a particle, while (Xj = when the site is empty. A particle configuration 
is collectively denoted by a = (o"j)igA. For each upward triangle j G A*, we denote the 
sites at the top, left and right of the upward triangle by t{j), and r(j), respectively. 
(See the left in figure [H) We describe a soft interaction among particles in j G A* by 
assuming a local energy Vj as follows: Vj{a) = when the particle number in j is zero 
or two, otherwise Vj{a) = J > 0. Explicitly, one may write 



Here, for a given site i, there are three upward triangles containing the site i. The 
upward triangles at the left, right and bottom of the site i are denoted by i{i), r{i), 
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and b{i) G A*, respectively. (See the center in figure ([T]).) We will see later that it is 
convenient to identify each upward triangle j with a dual site taken at the center of 
the same triangle. Then, the set A* may be identified with the dual lattice, which also 
forms a triangular lattice. See the right in figure [H 




Figure 1. Filled circles represent sites in the original triangular lattice A, and empty 
circles represent sites in the dual lattice A*. 

We study the equilibrium statistical mechanics with the Hamiltonian i/(cr) = 
J2jeA* ^■('^)- When studying statistical mechanics, it is convenient to express 

Hid) = ^Y1 [(2^*0) - - 1)(2^K.) - 1) + !]• (2) 

jGA* 

All the thermodynamic properties of the system with temperature T and density p are 
described by the Helmholtz free energy density 

/(T,p) = -T hm 1 log Z(T,p). (3) 

N^oo iV 

Here, the partition function is given by 

Z(T,p) = 5^e-^^('^)5(Arp,5^a,), (4) 

cr i 

where S{ , ) is Kronecker's delta function. Throughout the paper, 13 is identical to 
1/T, and the Boltzmann constant is set to unity. When we report results of numerical 
experiments, we always assume J = 1. We also define the ground-state entropy explicitly 
as follows. Let Q{E,p] N) be the number of particle configurations for a given density 
p and energy E. The entropy is defined by S{E, p; A^) = log fl{E, p; A^). By the ground- 
state entropy is meant Sgs{p', N) = min e S{E, p; N), and its density sgs(p) is defined 
as limjv^-oo >S'gs(p; A^)/A^- Obviously, sgs(p) > 0. Note that the free energy density 
/(T, p) is related to the entropy as /(T, p) = limAr_^oo(-E* — TS'(-E*, p; N))/N, where E^, 
is determined by dS{E, p)/dE\E=E^, = 

The spin model equivalent to the Hamiltonian ([2]) was introduced by Newman and 
Moore [13] and studied further [HI [15]. (The equivalence is exphcitly written in ([8]).) 
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These previous studies focused on the spin system without a magnetic field. As described 
later, this case corresponds to the lattice gas model with zero chemical potential or with 
half-filled density. A thermodynamic transition, which is our main concern, does not 
occur in this case. 

Ground states: The ground states of the model are defined as configurations in the 
limit T — > 0. Since H > 0, the configurations satisfying H{a) = are ground-states. 
Obviously, (Xj = for any i provides the ground state for the system with p = 0. As 
the other limiting case, the maximally packed ground states are easily found. They are 
crystals with p = 2/3. There are three configurations depending on its position. 

The other ground states whose density is neither nor 2/3 are given by non-trivial 
configurations. In order to see them, let us consider the case that one particle is put on 
a site for the vacuum p = 0. Then, the energy is non-zero in the three upward triangles 
containing this site. Thus, one particle must be put in each upward triangle so as to 
minimize the energy. This iterative process is not repetition, but it is given by a simple 
rule. 

We describe this rule explicitly. From the expression ([T]), one finds that the 
condition Vj{a) = is equivalent to the addition rule at(j) = cr^(j) + crr(j) (mod 2). 
Suppose that a particle configuration in the bottom line in the lattice is given. Then, 
when a particle configuration in the second bottom line is determined by applying the 
addition rule, the energy of upward triangles on the bottom line is zero. By repeating 
this process line by line, one obtains a particle configuration in A. This deterministic rule 
is equivalent to an elementary cellular automaton (CA) on a one-dimensional lattice. 
According to Wolfram's classification, it corresponds to Rule 102 included in class III 
[16], where space-time configurations generated by the rules in class III are known to 
be irregular for almost all initial conditions. 




Figure 2. (Color online) Space-time configuration of the cellular automaton. The 
blue circles represent sites at which at = 1. Periodic boundary conditions of the one- 
dimensional lattice (horizontal axis) are assumed, and time is the vertical axis. The 
initial conditions are given by a single seed at the center of the bottom line (left figure) 
and given randomly at the bottom line (right figure). 

As the simplest example, in the left of figure |2l we display a configuration starting 
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from just one particle in the bottom line. It is called a Sierpinski gasket which is a 
famous fractal pattern with the fractal dimension df = log 3/ log 2. Since the rule is 
linear, the superposition of two configurations (in the module two) also satisfies the CA 
rule. Thus, all the configurations generated by the CA rule are obtained by superposition 
of Sierpinski gaskets. In the right of figure [2], we display a configuration starting from 
particles randomly placed on the bottom line. 

Now, we consider the system with linear size L. That is, N = L"^. For simplicity, 
we first assume that the local energies for upward triangles whose bottom edge is at the 
top line or whose top site is at the bottom line are always zero. We also assume periodic 
boundary conditions at the side lines. With this choice of boundary conditions, all the 
particle configurations satisfying H{a) = are equivalent to space time configurations 
generated by L time steps evolution of the CA with the lattice size L. From the 
deterministic nature of the rule, the number of the configurations is 2^. Therefore, 
the density of ground-state entropy is estimated as sgs(p) < log2/L for any p, which 
leads to sgs(p) = in the thermodynamic limit. 

For the other boundary conditions, to determine ground state configurations for 
a given p is not easy. Concretely, let us consider the system with periodic boundary 
conditions in all directions. In principle, zero-energy configurations for a given system 
size L can be determined by an algebraic method [17]. However, for example, it has 
been known that the number of the configurations depends on system size L in a quite 
complicated manner which comes from the arithmetic nature of L. Furthermore, in some 
values of L, say L = 32, there is no zero-energy configuration when p ^ 0. In this case, 
the ground state configuration involves defects, each of which is defined on an upward 
triangle with the energy J. The best positions of the defects for minimizing the energy 
would be determined in an irregular manner as a function of p and L. Although the 
precise argument in this case is quite difficult, we provide a rough estimation of a typical 
number of defects as follows. We start with a simple observation that an isolated 
Sierpinski gasket with any size i is generated by assuming a defect at each vertex of 
the largest triangle. The number of particles involved in the gasket is l"^^. Let n be the 
highest level of the gasket included in the system. Then, 2" ~ 0{L). When the particle 
density is finite, the number of gaskets is 0(L^~'^f). Since there are three defects in each 
gasket, Ud — 0{L'^~'^^). Therefore, the ground-state energy is 0{L'^~'^^) and the density of 
ground-state entropy is estimated as sqs — 0(log2/L) +0{nd/N) ~ 0{L^^) + 0{L~'^^), 
which leads to sgs(p) = in the thermodynamic limit. 

3. Statistical mechanics 

Following a standard method in statistical mechanics, we consider the grand partition 
function 




(5) 



a 
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The pressure p(T, jj) in the model is given by 

P{T, /i) = T lim ^ logS(T, /i), (6) 
and the free energy density f{T,p) is obtained by the Legendre transformation 

f{T,p) = sup[/ip-p(T,/i)]. (7) 

Furthermore, by introducing a spin variable Si = 2o"j — 1 and parameters Ki = pj/2 
and K2 = -/3/i/2, we express S(T,/i) = e-^^'+^^'^^Y{Ki, K2) as 

s 

where s = (sj)jgA. The partition function Y is equivalent to that for the three-body 
Ising model on upward triangles under a magnetic field. Note that the case K2 = 
was studied in Refs. [131 [HI IIS] that the three-body Ising model on all triangles 
without a magnetic field was solved exactly [T8] . 

We start our analysis from the easiest case i^'i ^ 1 and K2 — 0(1), under which 
configurations minimizing J2jeA* M3)^Hj)^H3) contribute to Y. We assume boundary 
conditions such that all the spins outside the system are down. Then, the minimizer is 
uniquely determined as Sj = — 1 for all i. We thus obtain 

Y{K,,K2)c:^exp{{K, + K2)N) (9) 

as the leading order contribution. This result leads to p(T, jj) ~ 0. The first order 
correction to ([9]) arises from a single spin-flip with the probability exp(— 3/3J). This 
means 

p(T,/i) = 0(e-3/^^). (10) 

Duality relation: Next we derive a duality relation for p{T,p), which is presented in 
(l23|l . In the argument below, we assume K2 > 0. By making a change of variable 
Si —Si in the sum of (IE]), we have a relation 

Y{K,,K2) = Y{-K^,-K2). (11) 

With this relation, the expansion of the exponential function yields 



Y{K,,K2) = CJ2 



YI{1 + Si tanhK2) 



(12) 



JJ (1 + St{j)Se(j)Sr{j) tanhi^i) 
s IjeA* J \_ieA 

where C = (cosh i^i )^ (cosh ii'2)^- We then introduce a new variable aj whose value 
is either or 1 for j G A*. Concretely, dj = 1 when the term St{j)Si(^j)Sr{j) tanhKi is 
picked up in the expansion of the first product. Otherwise, aj = 0. In terms of this new 
variable, the first product is rewritten as 

JJ (1 + St(j)Se{j)Sr(j) tanhi^i) = JJ ist{j)S(,(^j)Sr{j) tanhi^i)'"^' 

jeA* jeA* 

= 5^(tanhKi)^.eA*^^ Jj4^«+^'-«+^^«, (13) 
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where a = {o'j)jeA* ^ind we have rearranged the product as that over i G A. Note that 
i{i), r{i), and b{i) were defined previously. See figure [H By substituting this expression 
into ( fT2l) . we obtain 



Y{Ki,K2) = C^^(tanhiri)^^6A*-. 



.iGA 



Here, we define by 

exp(-2ir*) = tanhfsTs, (15) 
and introduce a variable Sj = 2a j — 1. By an explicit calculation, we can confirm 

^^^^i(^)+^ri^)+^H^) + taufi Kg 4''"' ^^''''^ ) = 2e-^i* ^'''W ^^^0 + 1) . (16) 

Si 

Similarly, we write 

(tanhifi)^^eA*^^ = e-^2E,eA*(^-i+i)^ (17) 

where K2 has been defined by 

exp(-2is:*) = tanhfsTi. (18) 
By substituting (fT6!) and (fT7|) into ffT4|) . we obtain 

s 

where s = (sj)jgA* and A^* is the number of elements in A*. We have A^* = A^ when we 
assume periodic boundary conditions. By recalling that A* may be identified with the 
dual lattice, and by comparing (fT9l) with ([H]), one finds that the sum in ( |T9|) is nothing 
but the partition function of the same model with the new coupling constants and 
K2 defined on the dual lattice A*. See figured] again. The only and irrelevant difference 
is that the particles now interact with each other on downward triangles. We thus find 
the duality relation 

Y{K,, K2) = 2^(cosh A'i)^(cosh A'2)'^e-(^^+^2)A^y(i^*, k;). (20) 
We further simplify (fT5l), (fT8|), and (|20l) as 

Y{Ki,K2) = F(A:*, A:*)(sinh2A:lsinh2A:2)^^^ (21) 

and 

(sinh2A:*)(sinh2A:2) = (sinh 2A:*)(sinh 2A:i) = 1. (22) 

Finally, we define T* = 1//3* and /i* by K* = (3* J/2 and K* = -/3*fi*/2. Then, (|2T1) 
yields 

p{T, /i) = p{T*, ^*) - - + ^ + Tlog (2 cosh 4f cosh ) . (23) 



2T 2T^ 

Note that a duality relation connecting large (3 to small (3fi was discussed in Ref. [I9 
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Exact result for p = 0.5; As an application of the duality relation (1231) . we calculate 
the free energy density f{T,p = 0.5) exactly. When |yu| ^ 1, the dual point T* ^ 1 
satisfies exp(— J/T*) ~ —2p/T for any T. Since the pressure p{T*,ij;*) in the regime 
T* < 1 is known as we substitute it into Noting that 0{e-^^*-^) = 0(|/i|^), 

we derive 

J 



p(T,/i) = -^ + | + T 



log 2 cosh 



2T 



for \fj,\ <C 1. Since 



P 



dpiT,fi) 



we obtain from ([7]) 
f{T,p 



0.5) 



At=0 



J 



1 

2' 



-T 



log ( 2 cosh 



2T 



(24) 



(25) 



(26) 



This free energy density is equivalent to that of non-interacting Ising spins. 
Essentially the same result was obtained in Ref. [13]. Here, recall that the ground state 
of the system with p = 0.5 is given by superposition of the Sierpinski gaskets. Therefore, 
in contrast to the simple form of the free energy density, equilibrium configurations are 
complex as displayed in the left of figure [31 We then note that the equilibration time is 
extremely large, which might originate from complex configurations in the system with 
large /3. In the right of figure [3l we plot the heat capacity obtained numerically with 
the aid of an exchange Monte Carlo method [20]. The results for /3 > 4 still depend on 
the waiting time within our observations. This means that equilibrium states are not 
produced yet when f3 > 4. Indeed, the numerical data for /3 > 4 are different from the 
theoretical curve derived exactly. One may interpret that the system is in a dynamical 
glass phase. See Refs. [131 [HI [15] for detailed studies of dynamical properties in this 
regime of the model. With regard to this problem, we also remark that the mean-field 
model of the ferro-magnetic three-body Ising model without a magnetic field exhibits 
a glass transition in a meta-stable branch as well as a thermodynamic transition to a 
ferro-magnetic phase |21j . 



Thermodynamic transition: The most standard application of duality relations is to 
provide candidates of thermodynamically singular points. Indeed, if p{T, p) is singular 
at some point (To,/io), p{T,p) is singular at (Tq,Pq) too. Thus, a thermodynamic 
transition occurs on the self-dual line 

(sinh J/T) (sinh p/T) = - 1 (27) 

if there is one transition point in < T < oo for fixed /i < 0. Note that, in general, there 
is no reason to assume the latter condition, although there are many successful examples 
since the pioneering paper [22]. As one example, let us consider the case where T 3> 1. 
Then, the self-dual point is evaluated as /i ~ — TlogT. It is hardly expected that a 
transition occurs in such a high-temperature and largely negative chemical potential 
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Figure 3. (Color online) Numerical experiment of the model with ^ = 0. N = 1024. 
The left figure shows an equilibrium configuration for /3 = 3.75. The blue circles 
represent particles. The right graph is the heat capacity C as a function /3. The 
circle symbols represent experimental results, where the data for /3 > 4 still depend 
on a waiting time within our observation. The solid line is the exact result derived 
theoretically. 



regime. Thus, without a hasty conclusion, it is necessary to check whether or not 
transitions occur on the self-dual line. 

As another tractable case, we consider the system with < — /i/T ^ 1 fixed. For 
sufficiently small T, we have p{T,fi) ^ as derived in ( |TOl) . In this case, p ~ 0. Now, 
we increase T from the region p{T,fi) ^ (with keeping T <^ 1) and investigate the 
instability of the state p ~ 0. Following a standard phenomenological argument [23], we 
replace a macroscopic spherical domain by a spherical domain whose pressure is given 
by ( l24l) . The pressure difference from the bulk region is estimated by 



where we have used (I24p . If Ap > 0, the inserted domain grows if surface effects are 
ignored. That is, the first-order transition occurs at Ap = 0, which leads to 



The equality coincides with the leading order equation of fl27j) under an asymptotic 
condition that < —/i/T <C 1 and T -C 1. Therefore, we expect that the self-dual line 
under this condition provides actual transition points. 

In figure HI we display a result of numerical experiments for the system with 
fi = —0.04 fixed. It is seen that the density exhibits the discontinuous change from 
pi(p) to p2(p) at /3 = /3dis(p), where the left of figure H] shows that 2.75 < /3dis < 3.0. 
Since the dual relation (l27Il gives /3dis — 2.86 for p = —0.04, we judge that the theoretical 
consideration is consistent with the result of the numerical experiment. An equilibrium 
configuration in the lower temperature phase (T < Tdis(p)) is almost vacuum, while it is 
superposition of Sierpinski gaskets in the higher temperature phase (Tdis(p) < T ^ 1). 
Interestingly, as shown in figure [5], the intensity of density fiuctuations x = ^((p^ (p))^) 
measured with (T, p) fixed exhibits a power-law behavior x — (/3dis — in the higher 
temperature side, while no divergent behavior in the lower temperature side. Similar 





2T 



exp{J/T) 



1. 



(29) 
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Figure 4. (Color online) (Left): p versus /3 obtained by numerical experiments of 
the model with = -0.04. iV = 1024. (Right): The self-dual relation given by (P7)) . 
Note that ^ = 2.86 when fi = -0.04. 




Figure 5. (Color online) (Left): x versus /3 obtained by numerical experiments of the 
model with /i = —0.04. N = 1024. (Right): Log-log plot of x as a function of /3dis — /3 
for /3 < /3dis- The guide line represents x = Wdis ~ /3)~^^^- 

behavior was observed in several models [5]. Furthermore, we note that the transition 
temperature Tdis(/u) approaches zero as /i — )• 0. The equilibration becomes very hard 
as \fi\ is decreased. Indeed, we did not observe the first-order transition for the system 
with yu = —0.02 for N = 1024 within our observation time. In the left of figure M, we 
present a schematic phase diagram in the (/x, T) space, where < 0. Since there is no 
phase transition for sufficiently large T, the first-order transition line in the state space 
possesses a terminal point (/is, 71). 

Let us consider a schematic phase diagram in the state space (p, T). Since there 
exists a crystal configuration, a transition to a crystal state is expected to occur in 
the dense regime with p > 0.5. Indeed, we observed it by a preliminary numerical 
experiment. However, in this paper, we restrict our investigation to the regime p < 0.5 
in order to present a clear argument. 

We first recall that the density changes discontinuously from pi{p) to P2{p) at 
7dis(/^) in decreasing the temperature with /ig < < fixed. See the dotted fine 
in the right of figure [6l Then, the region below the two curves (pi(/i), Tdis(/i)) and 
{p2{p),Tdis{p)), with /is < /i < 0, is a coexistence phase that corresponds to the 
transition line in the left of figure [61 



Thermodynamic transition associated with irregularly ordered ground states 



11 



0.6 - 0.6 Disorder 




Figure 6. (Color online) Schematic phase diagrams in the state space (/i, T) (left) and 
in the state space {p,T) (right). In order to see the qualitative behavior more easily, 
the graphs are not accurately displayed. The chemical potential is equal to along the 
dotted line in the right figure. The coexistence phase is labeled by "Irregularly Order" , 
but its non-trivial nature is not ftrlly explored yet. See a discussion in the text. 

Now, we fix p satisfying < p < 1/2. There exists a special temperature Tc(p) 
below which the system is in the coexistence phase. Note that dp/dplr = holds in 
this phase. We further fix T to be smaller than Tc(p). The chemical potential p for a 
given (T, p) is determined uniquely. Then, from the analogy with a liquid-gas transition, 
one may imagine that two phases, (T, pi(p)) and (T, p2(p)), are separated. Note that 
pi ~ and p2 ~ 1/2 when T is sufficiently small(that is, p is sufficiently close to 0). 
Here, if this picture is correct, the interfacial energy proportional to 0{L) remains in 
the limit T — 0. However, as discussed in section [21 the ground states are far from 
such simple configurations. At the end of section |2l we conjectured that point defects 
appear with the energy 0(L^~'^f) in general ground states. Since this is less than the 
interfacial energy 0{L), we expect that phase separated configurations are not realized 
in the ground states. 

Of course, even if phase separated configurations do not appear in the ground 
states, this does not ensure the non-trivial nature of the system with finite temperature. 
We need to investigate the finite temperature system directly. More explicitly, we 
conjecture that the number of pure states in the coexistence phase is not two, but 
continuously distributed with a parameterization of density. That is, for example, when 
observing a density just on the first-order transition point in the (T, p) ensemble, its 
distribution function may form a continuous function, not a two-peaks function, in the 
thermodynamic limit. To prove this conjecture mathematically is a challenging problem, 
but careful numerical investigations will be done in future. With expecting the validity 
of our conjecture, we call the coexistence phase irregularly ordered phase. 

4. Future problems 

Before ending this paper, we list up problems to be studied in future. First of all, a 
precise phase diagram should be determined. In order to complete it, we will study the 
behavior near the critical point and the transition lines in detail. After the determination 
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of the phase boundaries, we will characterize each phase. Among the phases, the most 
non-trivial is the irregularly ordered phase. As an example, in figure U\ we present 
configurations for p = 0.25. The case f3 = 3.0 (left) is in the disordered phase, while the 
case /3 = 4.5 (right) is in the irregularly ordered phase. Both the configurations were 
obtained by long-time calculation of an exchange Monte Carlo method. The waiting time 
dependence is relatively small, but more careful checks are necessary to judge whether 
the system reaches the equilibrium state in the irregularly ordered phase. Numerical 
studies on this phase are quite challenging. 




Figure 7. (Color online) Configurations for the system with p = 0.25 fixed. The blue 
circles represent particles. /3 = 3.0 (Left) and /3 — 4.5 (Right) 

Theoretically, the next step study is to construct a perturbation theory for the 
system with small T and small |/i| in order to understand peculiar features observed 
numerically. This study will formulate calculation techniques for treating excitations 
of fractal patterns. Here let us recall that the mapping to a free gas of defects was 
useful to understand behaviors for the system with fi = [HI [15] . It is then a natural 
idea to describe the small system through a week interaction among defects. For 
example, it might be interesting if one confirms that defects form a bounded state in the 
irregularly ordered state. Furthermore, the study on algebraic properties of the cellular 
automaton rule may provide a hint for characterizing the irregularly ordered phase. As 
a different direction of study, one may investigate the irregularly ordered phase from a 
viewpoint of thermodynamic glass. At present, the existence of thermodynamic glass 
in finite dimensions is not established yet, mainly because dynamical singularities cover 
the equilibrium nature. Nevertheless, the understanding at the mean-field level has been 
greatly accumulated. In particular, there is consensus that thermodynamic glasses at 
the mean-field level are characterized by a one-step replica-symmetry breaking [2H |25] . 
Recently, a thermodynamic glass phase in a lattice model [26] has been investigated 
numerically along with this idea [27]. As did in the study, one might argue a replica- 
symmetry breaking in the irregularly ordered phase of our model. 

Finally, we remark that our model has no direct correspondence with experimental 
systems. However, needless to say, it is quite amazing to find phase transitions associated 
with irregularly ordered ground states in natural phenomena. (See Ref. [28j for some 
connection of the model to experimental systems.) We expect that the phenomenon 
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reported in this paper will be studied further theoretically and experimentally. 
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